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Abstract 

We construct non dispersive two soliton solutions to the three dimensional gravita- 
tional Hartree equation whose trajectories asymptotically reproduce the nontrapped 
dynamics of the gravitational two body problem. 

1 Introduction 

1.1 Setting of the problem 

We consider in this paper the three dimensional gravitational Hartree equation: 
iu t + Au — 4>\ U \2U = 0, 

A<p lul 2 = \u\ 2 i.e. </> H 2 = -j^pij * \u\ 2 , (1) 
(t, x) G R x R 3 , u{0, x) = u Q (x), u : R 3 -> C. 

This system arises in Physics as an effective evolution equation in the mean field limit of 
many body quantum systems, see for example [SJ, [7j. 

It is well-known, see [3] and references therein, that the Cauchy problem for (pQ) is 
globally well-posed in the energy space H 1 = {u,Vu G L 2 (R 3 )}, i.e., for u G H 1 , there 
exists a unique global solution u(t) G CiR,!! 1 ) of ©. 

Moreover, the following quantities are conserved by the H 1 flow: 



I? norm: J \u{t,x)\ 2 dx = J \uo(x)\ 2 dx, 

Hamiltonian : H(u(t,x)) = — J \X7u(t,x)\ 2 dx — — J |V0| u |2(t, x)\ 2 dx = H(uq) 
Momentum : Im (^j Vuu(t, x)dx^j = Im (^j X7uoUo(x)dx^J . 
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Recall also that the conservation of the energy and the 1? norm imply an a priori bound 
on the kinetic energy as derived from the Hardy-Littlewood-Sobolev inequality and then 
the Gagliardo-Nirenberg inequality: 

Vu G H\ J |V0 H 2| 2 < C||w||*^ < C\\Vu\\ L2 \\u\\l 2 . 

This implies that (H|) is subcritical so that all H 1 solutions to (H|) are global and bounded 
in H 1 . In contrast, note that in dimension four, the Hartree equation with potential jj-^jp 
is I? critical and a singularity formation is known to possibly occur, see |13| for recent 
progress on this problem. 

Equation ([1]) possesses a large group of symmetries: if u(t, x) is a solution to (P), then 
for all (t , a Q , fa, 70, A ) € R x R N x R N x R x R+, so is 

v(t,x) = \ 2 u(\ 2 Q t + t ,\ x + a - foty^-^e^ . (2) 

Special solutions are expected to play a fundamental role for the description of the 
dynamics of ([I]). They are the so-called solitary waves, of the form u(t,x) = e lt W(x) 
where W(x) solves 

AW - (j)\ w \2W = W. (3) 

We shall denote by Q the so-called ground state solution to (j3J), which is defined as the 
unique radially symmetric nonnegative solution to (|3J). The existence and uniqueness of 
the ground state has been proved using variational and ODE techniques by Lieb [16J. 

The group of symmetries ([2]) then generates an eight parameter family of ground state 
solitary waves: 

Qa,p,y,\(t, x) = \ 2 Q(\x + a - ptyk^-l ')e iA2 *e% (a, p, A, 7) G R 3 x R 3 x R x R* + . 

Observe that the center of mass of such a solution is evolving according to the free Galilean 
motion with constant speed (3. Recall that a variational characterization of the ground 
state Q has been derived by Lions [17] using the concentration-compactness technique. 
Given M > 0, the minimization problem 



inf H(u) 

II"" 

is attained at 



M L 2=M 



u{x) = \ 2 {M)Q{\(M)x + a)e i7 , (a, 7) G R x R, \{M) = 



2 



Moreover, every minimizing sequence is relatively compact in the energy space up to phase 
and translation shifts. Following Cazenave and Lions [5] and using Lieb's uniqueness result 
|16j . this automatically implies the orbital stability of the ground state solitary wave. 
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The question of the long time dynamics of the Hartree equation is widely open. How- 
ever, in other settings like the one dimensional (gKdV) equation or the Schrodinger equa- 
tion (NLS) with power nonlinearity 

f iu t = - Au - M p-1 u, (i.ijelxl" 

{ u(0,x) =u (x), u :R N ^C, p<l + j r , U 

multisolitary wave solutions are conjectured to be the building blocks for the description 
of the long time dynamics. Roughly speaking, as time goes to +00, a generic solution 
should split asymptotically into a sum of solitary waves which move away from each other 
and a radiative term which disperses. There is a very extensive literature on this subject, 
we refer for example to [3], [25], [20], [21], [6], [28], [27] (and references therein). 

It seems that one of the first problem to address for the understanding of the long time 
dynamics is the existence of non dispersive multisolitary waves. Apart from completely 
integrable systems like the (KdV) equation or the cubic one dimensional (NLS) where 
explicit multisolitary waves can be exhibited (see e.g. pD]), the existence of such objects 
goes back to Merle [23] for critical (NLS) problems, and has later been extended by Martel 
|18j to subcritical (KdV) equations and by Martel and Merle [19] for general subcritical 
nonlinear (NLS) problems. See also Buslaev and Perelman [3] and Rodnianski, Soffer and 
Schlag [25] for different approaches. 

Let us observe that for a system like the subcritical (NLS) (J2J), the sum of two solitary 
waves moving away from each other at constant speed is a solution to ((4|) up to an ex- 
ponentially small in time correction. This suggests that there is essentially no interaction 
between the two ground states, and that nondispersive multisolitary wave solutions can 
be constructed with asympotic trajectories given by the free Galilean motion along any 
chosen non parallel lines, [19] . 

1.2 Statement of the result 

This paper concerns the construction of two soliton solutions of the Hartree equation. Note 
that for the Hartree equation, in contrast with the (NLS) or (gKdV) case, even though 
the ground state is exponentially decreasing in space, the nonlinearity is long range from 
the slow decay of the gravitational field 

j. t ^ C 

^IQI 2 ( r ) ~ — as r +°°- 

The consequence is a strong coupling between the two solitons. In this paper, we claim 
that this constraint implies a non trivial dynamics for the centers of mass of the solitons, 
which - as one should expect from physical grounds - asymptotically converges to the 
dynamical system of the two body problem in Newtonian gravity. 
Now, we recall the well-known definition of the two body problem: 
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Definition 1 (Dynamics of the two body problem) A trajectory of the two body 
problem with masses (Ai,A2) £ (M+) 2 is a solution (a\(t), 012 (t), /3\(t), fait)) to 

(5) 

P = fo -Pi- 

The center of mass evolves according to the free Galilean motion, i.e. \2ct1 + \\dt2 = 
and a(t) evolves in a fixed plane. Moreover, the dynamical system admits the conserved 
Hamiltonian: 

\Ql\h / 1 . 1 \ 1 




The dynamics are classified as follows: 

(i) Hyperbolic trajectory: If Eq > then a(t) describes a hyperbola with 



lim II^lZh = 2^0- (6) 

t— >+oo t 

(ii) Parabolic trajectory: If Eq = then a(t) describes a parabola with 

lim = C > 0. (7) 

(Hi) Elliptic trajectory: If Eq < then the dynamic is periodic in time and a(t) 
describes an ellipse. 

In the rest of this paper, we fix (Af^A^) S (R+) 2 , and we consider a given solution 
(af (i), af (t), f3f°(t), 0f{t)) of the two body problem with masses (Af , \f). We assume 
throughout the paper that the center of mass is fixed at the origin and that the trajectory 
lies on the plane {#(3) = 0} where x = (xm, %(2)> x (3)), 1 - e - 

Vt > 0, A^a^(t) + Af a^i) = and (a°°) (3) (t) = 0, (8) 

where aj(t) = ((oij)m(t), (ay)(2)(i), (ay)(3)(£))- Note that these assumptions do not restrict 
generality by standard use of translation, rotation and Galilean invariances. 

Our main result concerns the existence of non dispersive two soliton solutions which 
asymptotically reproduce the non trapped dynamics of the two body problem in the fol- 
lowing two cases: 

1. Hyperbolic case without restriction on the masses; 

2. Parabolic case with equal masses. 

Theorem 1 (Existence of a two soliton with a non trapped trajectory) 

Let (Af.Af) G (R* + ) 2 and let {af{t),oif{t),j3f{t),l3f{t)) be a solution to the two body 
problem with masses (Af^A-f) such that |2P holds. 
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1. Hyperbolic case. Assume that (af 3 ^), a^(t), P%°(t), /?2°(*)) is hyperbolic (Eq > 0). 
Then, there exists an H 1 solution of (CP and 7i(£),72(i) such that 



lim 



u(t, x) 



2 



1 



-Q 



af(ty 



0) 



H 1 



2. Parabolic case. Assume that (af°(i), a^(t), f3f°(t), 02°(t)) is parabolic (Eq = 0) and 
= = A°°. Then, there exist an H 1 solution of |iy and ^(t), 7i(*)> 

72 (t) such that 



lim 



3=1 



1 n f ^^iWA ^(Wtth 



0. 



H 1 



lim > 



t— >+oo 



fe=l 



(«i)(fc)(*) 



+ l(«i)(3)(t)l = J i = 1,2. 



Comments on the result: 



1. Constraints on the asymptotic trajectory. Note that for local nonlinearities with 
weak interactions, like in the case of power nonlinearity (NLS) problems, the asymptotic 
dynamics of the centers of mass can be prescribed along any prescribed lines, see |19j . In 
contrast, in Theorem[H the asymptotic hyperbolic trajectory has to be along two coplanar 
lines. We expect this situation to be the only possibility. 

2. Parabolic case. The behavior displayed by the solution constructed in Theorem 1 in 
the parabolic case does not correspond asymptotically to an almost free Galilean motion 
of each body, i.e. along two coplanar lines, which is for one soliton the natural motion 
from the Galilean symmetry. Thus, the result of Theorem 1 in the parabolic case displays 
a new nonlinear (critical) regime where the trajectories of the two solitary waves do not 
correspond to asymptotically free solitary waves. 

Note also from the proof that we do not know whether ||ay(t) — otj°(t) \\ — > as t — > +oo 
in the parabolic case. This question is related to the computation of a constant, which 
does not seem to be explicit. However, if we set V°° = {a < j )0 (i), a^°(i), t > 0}, then 



lim dist((a 1 (t),a 2 (t)),V c 

t— >+oo 



0, 



which means that (ai(i), ct2(t))) asymptotically describes the same set as (af 3 (t) , a?? (t)) 
with different speed. See Appendix A. 2. 

3. Connection to existing literature. The fact that the Hartree problem displays 
solutions which reproduce in some sense the point particle newtonian interaction is not a 
surprise. This has been in particular observed by Frohlich, Yau, Tsai [9] where effective 
modulation equations leading to the Newtonian Galilean motion for the center of mass 
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are derived in certain asymptotic regimes - but not for all times t — > +00. In [TO], Gang 
and Sigal describe the long time behavior of the (NLS) equations in a trapping very well 
localized potential and show how the Newtonian law of motion is here again driving the 
center of mass of the ground state part of the solution. 

4- Uniqueness. For the hyperbolic case, one may ask whether a solution satisfying 
is unique. By the method of this paper, one can expect a weaker statement, i.e. uniqueness 
in a class smaller than ([9]). Uniqueness in general is an open problem, certainly related to 
the -delicate- stability problem, see next comment. See |18| for a general uniqueness proof 
in the (gKdV) case. 

5. Stability. The main open problem after this work is the stability of the two soliton 
solution as derived for some (NLS) and (gKdV) type problems, [25], [20], [21]. We expect 
the instability of the parabolic trajectory and the stability of the hyperbolic one. Another 
important open problem concerns the non existence of trapped elliptic type trajectories 
which most likely are destroyed by the dispersive effects of the flow. 

The proof of Theorem Q] relies on a refinement of the techniques developed by Martel 
and Merle [18], [TO], and proceeds into two steps. 

• Construction of an approximate solution to all orders: The first step is to display 
a general framework to compute an approximate solution to the two soliton equa- 
tion with an arbitrarily high order error jj^jpv where ||a|| is the distance between 
the center of masses. Here our analysis relies on a finite dimensional reduction of 
the problem. This kind of procedure is reminiscent to the computation of formal 
ansatz and solutions in various nonlinear dispersive settings including for example 
the description of singularity formations, [26], [23], [14] , the description of the long 
time dynamics of the (NLS) near a ground state, [TO], or the interaction of two 
solitons [22]. We propose here a systematic way of producing such approximate so- 
lutions. The main advantage is that the approximate solution is now good enough 
to somehow reduce the problem to an almost short range problem. 

• Construction of the nondispersive two soliton: The second step is to build the exact 
two soliton solution by solving the problem for the small remainder backwards from 
infinity with a procedure which is reminiscent from the construction of nonlinear 
wave operators -see for example [2] for a similar strategy-. The control of the re- 
mainder relies on a space localization of the conservation laws and energy estimates 
on the linearized flow close to a solitary wave. Such estimates are a consequence 
of the variational characterization of Q as first exhibited by Weinstein [30] in the 
local (NLS) setting, and some elliptic non degeneracy properties proved by Lenzman 
|15] for the Hartree problem. While this strategy has the advantage of not requiring 
any dispersive estimate, the localization of conservation laws creates large errors in 
our setting due to the long range structure of the problem and the only polynomial 
distance between the two solitary waves. In fact, it turns out that the closure of 
these estimates is critical here with respect to the £ decay of the gravitational field, 
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in particular in the parabolic case where we can only treat the case of an asymptotic 
symmetric two body problem. Eventually, the success of our strategy relies in a cru- 
cial way onto the the construction of an arbitrary large order approximate solution 
from the first step, see Lemma [7] and Remark [U 

We expect this strategy to be quite robust and to apply to a very large class of situations 
and problems. 

Acknowledgements. J. K. is supported by NSF-Grant DMS-757278 and a Sloan Fel- 
lowship. Y.M is supported by ANR Projet Blanc OndeNonLin. P.R is supported by ANR 
jeunes chercheurs SWAP. Part of this work was done while J.K was visiting the Institut 
de Mathematiques de l'Universite Paul Sabatier, Toulouse, which he would like to thank 
for its kind hospitality. 

2 Construction of the approximate solution 

The purpose of this section is to construct a formal approximate two soliton solution of 
([T]) at any order of tt^tt, see Proposition [T] and Proposition [2j Our strategy relies on a finite 
dimensional reduction of the problem which is connected to the strategies developed for 
example in [26], [21], [10], [22], [H]. 

2.1 Derivation of the coupled equations 

Let us start with computing the evolution equation in renormalized variables. Let 




Then 



idtu + An — (f)\ u \2U = —j i\ 2 dtv + Av — v — (f>\ v \av 
A L 

-iXXAv - \ 3 {P ■ y)v - iX(a - 2/5) • Vv 



(10) 




where A denotes the differential operator 



Av = 2v + y ■ Vv. 



We compute an ansatz of a two soliton solution as t — > +oo by letting 



u(t, x) 



ui(t,x) + u 2 (t,x) 



+ 



i 



v 2 (t,y 2 )e- i ^e lf3 ^- 
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where 



The nonlinear term is 



x — a jit) 



■\u\ 2U = ^|Mi+n 2 | 2 («l + u 2)- 



(11) 



By construction, u\ and u 2 will have disjoint supports up to exponentially small correc- 
tions. Hence we will not consider the crossed term for now (see Proposition [2] for details). 
By rescaling, 

3 



and so 



4>\ U \2U = 0| Ml p + 0| U2 |2)ui + (0| Ul |2 + </>|„ 2 |2)u 2 + 0(e 7t ) 



A? 



1 

+ 53 



«1 



A 2 ) 

2 



|2 + I — I 0|„ 2 |2 ( —2/1 



Ai a 2 — oi 
A 2 



Ai7 



A 2 

A 2 a 2 — qi 



A, 



Thus, we rewrite the equation 



ad t u + An - H2 u = ^ rj^^aJe-^We^W' 1 



with 



Eqj(t,x) = iXjdtVj + Avj 



a 2 — ai 



Aj+i Aj + i 



iXjXjAvj — A|(/3j • — iXj(dj — 2(3j) ■ Vvj 



and the convention of mod(2) summation: 

j + 1 = 1 for j = 2. 



2.2 Dipole expansion of the gravitational field 

Let 

a(t) = a 2 (t) - ai(t), /3{t) = /3 2 (i) - ft(t). 



(12) 
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We expand the gravitational field created by the soliton Uj+i on uy. 



i+i 



i+i 



A 



•i+i 



4vrA 



\a- (-iy{\ j+1 £- XjVj] 



We proceed to the formal dipole expansion of the gravitational potential which main 
contribution corresponds to the region \\yj\\ + ||£|| <C ||a||. Such a formal argument will be 
justified in the proof of Proposition [2] below for functions Vj with sufficiently fast decay at 
infinity in space. 
Let 

1 1 



ll«-CII 

and use a Taylor expansion: 



a 



C-at IICII 2 



\a 



\a 



(13) 



\tt 



1 - / 



CI 



N+l 



O 



\N+1 



(14) 



Explicit computations yield for the first terms: 
1 



F 1 (a,C) 



\a 



i? 2(a,C) = -^|, F 3 (a,() 



a 



a 



a 



a 



Observe in particular that is homogeneous of degree —k in a. 
We thus introduce the approximated correction field of order N: 

N 

k=i 



(15) 



with 



iN,k,app ( n 



|r; i+1 (0| 2 ^((-l) j «, A,+ie - A iyi )de- 



(16) 



2.3 Finite dimensional reduction of the dynamics 

From the previous formal calculations, we aim at constructing an approximate solution to 



the system: Eq 



(AT) 



for j = 1, 2, where 



(AO 



iXjdfVj + A 



...12+01 |2 



iXjXjAvj — ■ Uj)vj — iXj(ctj — 2/3j) ■ Vvj 
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We now proceed to a finite dimensional reduction of this problem and look for a solution 
of the form 



AN) 



(AO 



(a(t),/3(t),Ai(t),A 2 (t)) 



Here vj N ^ is stationary and hence Vj N ^ depends on time only through the modulation 
parameters t — > (a(t), /3(t), Ai (t), \2(t)) to be chosen: 



d t v 



(JV) 



E 

fc=i 



(AO 



<9A fc 



-A fc + 



(JV) 



da 



a + 



We obtain 



AV' (A '» - V m 



iN,app 
W\2 + V^AN)^ 



W"'\ 



i\jM^Avj N ^ - \pf< ■ y jV] 



(JV) 



V 



(JV) 



2 oy(AJ) 



(JV) 



fe=l 



<9A fc 



da 



BV (N) 



17) 



<)) 



(JV) 



where S 1 -^ encodes the finite dimensional system of the geometrical parameters: 



- - - l X j ia j -2p j ).VV^-X^ J -B^).y ] vW 



- iX 3 (X, - M^)AVP + A|(7, + ^ - H/3,11 2 - 4- • ay)*f> 

(iV) ay (AT) 



k=l 



(X k -M { k N Y-^— + (-l) k (a k -2(3 k ) 



dV 



dX 



da 



(AT), ^ 



(JV) 



9/3 



(18) 



Now, we build an approximate solution using a series expansion of vj N ^ in (a, (3, Xi, X2) 
with a prescribed structure. We need the following definition. 

Definition 2 (Admissible functions) Let n € N. 

(%) H^e define S n the class of complex valued functions a(a, (3, Ai, A2) : (K 3 \ {0}) x (]R 3 \ 
{0}) xR* + xR* + ^C of the form: 



n3 



a{a, (3, Ai, A 2 ) = X] x i X 2 ^(2)^) fhMMMM^)^ "1,^2, n 3 € N, 



2 Ji,i 2 ,fcl,fc2,fc3=0 
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where (3 = (Pm, P(2)> P(3)) an d fhMMMM ^ s a complex valued C°° function of a away 
from a = and homogeneous of degree —n in a. Then, we set 

deg(a) = n. 

(Hi) We say that a function tp : (a, (3, Ai, X 2 ,x) : (l 3 \{0})x(I 3 \{0})xi;xR;xf 3 
is admissible of degree n if it admits a finite expansion of the form 

"4 

rp(a,/3,X 1 ,X 2 ,x) = ^2 o- m (a,(3,X 1 ,X 2 )T m (x), ra 4 G N 

m=0 

where a m G S n and T m is well-localized in x: 

Vk > 0,Vx G M 3 , V k T m {x) < e ~ Cm - kM . (19) 

Then, we set 

deg{ip) = n. 

(iv) We say that a function ip : (a, j3, X u X 2 , x) : (M 3 \{0}) x (8 3 \{0})x1;xK;xR 3 -> C 
is admissible of degree > N if it admits a finite expansion of the form 

1p = Ipn, n 5 > N 
n=N 

with ip n admissible of degree n. We say that deg{ip) > N with abuse of notation. 
We claim the following technical facts. 

Lemma 1 Let {n\,n 2 ) G N 2 . Letip\,ip 2 be admissible of degree respectively n\,n 2 . Then: 

(i) ip\ip 2 is admissible of degree n\ + n 2 . 

(ii) fy^ipi is admissible of degree n\ + n 2 . 

(Hi) ViV > 1, p^ app ip 2 is admissible of degree > 1 + n\ + n 2 . 
(iv) VI < k < N, 4>^ k,app ip 2 is admissible of degree k + n\ + n 2 . 

Proof. Properties (i)-(ii) are clear by the definition of admissible function of degree 
n (see Definition [2]) and the properties of homogeneous functions. Property (iii) is a 
consequence of (iv). 

Finally the proof of (iv) relies onto the fact that F^((— l) J a, Aj + i£ — Xjyj) appearing 
in (|16h is a polynomial in (Ai,A2) which coefficients are homogeneous of order —k in a 
and k in £. Indeed, -Ffc(a, C) is the k-th term to the Taylor expansion (|14p of (|13p . 

Then, we claim the following Proposition. 
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Proposition 1 (Expansion of the approximate solution) Let Ai,A2 > 0, {a, pi) G 
(M 3 /{0}) 2 and N > 1. Let j = 1,2 with the convention j + 1 = 1 for j = 2. We can find 
an asymptotic expansion 



v} N \ yj ) = Q{yj) + £ 2f >(a, /?, Wl Al) A 2 



iV 



(20) 



n=l 



TV TV 

M\ N) (a, /3, Ai , A 2 ) = £ mf> (a,/3,X 1 ,X 2 ), B^ N) (a, (3,Xi,X 2 ) = J2 bf ] (a, p,\i,X 2 ), 

n=2 n=2 



such that the following holds true. 

(i) Tj is admissible of degree n. 

(ii) mf eS n , bf £S n . 

(ii) Approximate solution of order N + 1 : Let 



(21) 



3 3 



i . ,N,app 



V, 



(TV) 



(22) 



+ i 



' 2 8V {N) BV {N) 



k=l 



dX k 



da 



2(3 + 



dV 



(TV) 



dp 



Then Eq^ is admissible of degree > N + 1. 

The proof of Proposition [1] proceeds by induction on N . We start with computing 
explicitely the first two terms in powers of ||a|| . 



2.4 Computation of the first terms 

First, we recall the following standard facts on the linearized operator around Q 

L+f = -Af + f + 4> Q2 f + 2cf> {Qf) Q, L_/ = -Af + f + <t> Q 2f. 
which are consequences of the nondegeneracy of the kernel of L + , L_ as exhibited in |15| . 

Lemma 2 (Inversion of L±) Let n G N and f be real-valued, admissible of degree n. 

(i) If (/, VQ) = ; then there exists a solution to L + u = f which is real-valued and ad- 
missible of degree n. Moreover, if f is radially symmetric then u can be chosen radially 
symmetric. 

(ii) If (/, Q) = 0, then there exists a solution to L_u = f which is real-valued and admis- 
sible of degree n. Moreover, if f is radially symmetric then so is u. 
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Proof, (i) Introduce the Hilbert space TL of all real- valued functions / G L 2 (M 3 ) which 
satisfy the orthogonality relation 

(/, VQ) = 0. 

Then note that L + maps TL D D into TL, D being the domain of L + in L 2 (IR 3 ). It is 
proved in |15j that the kernel of L + consist precisely of the components of VQ. Hence for 
v G TL PI D, we have L + v = iff v = 0. Now consider the equation 

L+u = f 

for given / G TL. Writing 

we can write the above equation as the following problem: 

v + Kv = f, K = K 1 o(-A + I)-\ K lV = [(f> Q 2V + 2(f> {Qv) Q}. 

The operator K : TL — > TL and is compact. Thus, its adjoint K* : TL — > TL, is also compact. 
Moreover, from the information on the kernel of L + , I + K* has a trivial kernel: indeed, 
denoting the projection in L 2 (M?) onto (span(d xl Q , d X2 Q , d X3 Q)) ± by II, we have 

k* = no (-A + rr 1 o/^ 

Hence the assertion (/ + K*)v = for some v £ TL implies 

v + ((-A + I)- 1 oK 1 )v = a-VQ 

for suitable a£l 3 , Standard elliptic theory implies v G C°°(1R 3 ), and then 

L+u = (-A + /)(«• VQ). 

Taking the inner product with a ■ VQ yields 

= ((-A + I) (a ■ VQ),a • VQ), 

which immediately implies a = 0, whence L + v = 0. But since v G 7Y, the result in [TS] 
implies v = 0. Using now the Predholm alternative for / + K* , we see that this operator 
is invertible. Thus I + K is also invertible. 

In order to obtain exponential decay for u, we use the well-known argument of Agmon 
PQ. For f & TL and < c < 1 assume a bound < e~ c " a: ". Write the equation above 

with u G TL as 

u = (-A + I)" 1 / - (-A + iy 1 o K lU 

Note that standard elliptic bootstrap methods imply u G -ff 2 (R 3 ) and in particular u G 
L°°(]R 3 ). If we use the explicit kernel representation for (A + I) , we get the equation 

r e -\\ x -y\\ r e -\\ x -y\\ 

7rs 4vr||x - y\\ J R3 4tt\\x - y\\ v W) 
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Now choose < 5 < \c such that Q(x) < e s » x \ and pick a large M >> 1. We now 
estimate minjilf, e* 5 '^" }u(x). First, we note that 



\\ x -y\\ 



4:ir\\x — y 



-J{v)dy 



< C. 



Next, choose L large enough such that 



-(i-<5)INI 



: E l<My)l)ll ^R 



< 



Then we get 

|min{M, e 6M } 



2U + 2cj> Qu Q}(y)dy\ < \ sup[min{M,e«}Ky)|] 



A-k\\x — y\ 



+e 



S\\x\\ 



\\ x -y\\ 



\ y \<L 4tt\\x - y\ 



-6 Q 2udy\ + 



IF-Vll 



A-k\\x — y 



;2<t>(Q U )Q(y)dy\ 



1 



< - aup[mm{M,e SM }\u(y)\] + C, 

2 , 



where C is independent of M. Summarizing, we have shown that 

min{M,e 5||a;|l }n(x) < C + - sup [min{M, e SM }\u(y)\] 



with Co independent of M. Since sup J/gK 3[min{M,e 5 " 3/ "}|u(y)|] < oo, taking the supremum 
over x € M 3 yields 

sup[min{M,e 5 " a:|l }|u(a;)|] < 2C , 

xeM. 3 

with Co independent of M. Letting M — > oo, we obtain the exponential decay of u(x). The 
corresponding statement for the derivatives follows analogously, provided the derivatives 
of / decay exponentially. 

The statement concerning the admissibility is obtained by linearity, and (ii) is proved 
like (i). Parity preservation is also is standard fact. This concludes the proof of Lemma [2 

Now, we proceed to the computation of the first terms of the approximate solution. 
This will provide an initialization of the induction argument of the proof of Proposition [TJ 
Moreover, the exact expression of the first two terms is required in the rest of the paper. 



j2,app 



\ 2 

47rAo_i_i \\a 



-i-iy 



\ 2 



47rAj + i \\a\ 



Aj+i 



(a -S + i(0l 2 *-A>' a) 



a 



\v j+ x(0\ 2 dt; 
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Order JL: We choose mf = bf> = and vj 1] = Q + T^ 1] such that: 

-L+lf + J^L Q = 0. 



47rAj + i ||a| 



Note that (Q, VQ) = and thus the existence of an admissible, real-valued and even 
function 2^. of degree —1 follows from Lemma [2j 

Order jrKn'- We set V- 2 ^ = Q + T- + Tp^ where T- 2 ^ is complex- valued and matches 

IHI ^ 3 3 3 v 

the O(-p-jp) terms. First observe from the radial symmetry of Q that J Q 2 (C)(C ' ce)d£, = 
and thus the order 2 term in the dipole expansion is exactly 



HQ" 2 " l j ^4vrA J+1 ||a||3 J l ° 



ST? 1 ' (2) (2) 

through the term — • 2/3. We thus look for m) , b) such that we can solve the 



Moreover, T- generates a contribution of degree (—2) to the imaginary part of Eqj 

through the term 
following system: 

- L_(/m(T) z 0) + A j^~ • 2 /3 + A/m^AQ = 0, (23) 



,(2)» | ^2 C ' i j no , x (2), 



-L + (i?e(i; (2) )) - ^dtW - <l> {T m )2 Q 



47rAj + i||a 



||Q||| 2 t£\ + 2(7 QT^Q 



J'^IL l(«-yi)Q = o- (24) 

From Lemma O we may solve (I23p with Im(Tj ) admissible of order (—2) provided we 

(2) 

adjust the scaling parameter such that 

- A,mf (AQ, Q) = -Aj- J!^mf = A 2 ■ 2(3, q\ (25) 



so that rrij ' G 5 2 



(2) 

(2) 



Similarily, we compute b, in order to have that the inner product of (|24j) with VQ is 
zero, and this yields: 

2 

3 v ' 47rAj + i 1 1 ck | 



h m = { _ iy+l \^a 2 _ (2g) 
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(2) (2) 

Note b 2 — b\ is the first term in the formal expansion of j3. Thus, (|26p . written for 



j = 1,2 corresponds to the modulation equations on the relative position a 

= 02 -Pi' 

1 I \ n 

+ 



OL2 — Oil, 



2/3, P 



WQWb 



\a\ 



4ir \Xi \2 / 

that is the two body problem interaction in Newtonian gravity with masses t~, tt- 



Remark 1 From 



and 



, we infer that Re(Tj 2 ^) is even, since the term in 

(2K 



completely eliminates the only odd term in the equation of Re{T, ). 



Remark 2 Let us observe that there is no uniqueness in the choice ofTj'^ as one can 
always add non trivial elements of the non trivial kernels of L + ,L-. However, this would 
in turn modify the modulation equations and then one can explicitely check that in fact, the 
same full function Vj(t,x) would be produced in the end. Hence an asymptotic expansion 
of this form to high order is indeed unique. 



2.5 The main induction argument 

We now proceed through the proof of Proposition Q] by induction on N. 

Proof of Proposition Q] 
Step 1. Expansion of the order N+l error. 
We argue by induction. The case N = 2 has been completely described in the previous 
subsection. We assume ./V and prove N + l. Let V, 



M^ N) + m 



(N+l) B [N+l) 



(N+l) _ y (N) T (N+1) M (N+1) 
3 J ' 3 



Bj N ^ + bj N+1 ^ . Then an explicit computation gives: 



Eq 



(N+l) 



L + (Re{Tf + ^)) + ^(b {N+l) 



3 " 1 ■ J\ 3 



L^LmiTf^'^ + Xjim)' 



) ■ yjQ 

)AQ 



+ Eq 



(N) 



+ Err^p + iErrp + Errp + Errp 



Errp encodes linear type errors on the geometrical paramaters: 



;JV+l,JV+l,OPP 

1QI 2 



Q 



Err 



(i) 



-i\i 



M 



(N) AT (N+1) + {m (N +l) )A ^ T j 



N+l 



3 



k=l 

V3 E T ^ k) 

k=l 



3 



(27) 
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(2) 

Err- is the error due to the formal time derivation: 



ErVj 



(2) 



^ 2 E 

k=l 
Q V (N+1) 

dx 



(N+l) 



f)T [ 



(N+l) 



(b 



dXk 



(N+l) ,(N+1) 



dT, 



(N+l) 



dT, 



(N+l) 



+ 



+ 



dp 



(B. 



(N) 



da 



2(5 



(28) 



(3) 

Err - encodes the nonlinear type errors in the expansion of the gravitational field created 
by Vj : 



Err 



(3) 



,(JV+1) 



+ 2(j> 



Re(vj N) ,T^ N+1 >V' 



+ 



Re((v} N) -Q)X 



(N+l) 



T),Q + <l>\ T (N+l)\2 V i 



(N+l) 



(N+l) 



Q) 



I 2 J 



(29) 



Err^f 1 encodes the nonlinear type errors generated by the gravitational field cjh ,app : 



Err {4) 



iN,app 
| T (JV+1)| 

\ 1 j+i 



+ <P 



N,app 



2Re(V^>,T, 



V: 

(JV) ^(iV+lK I 3 



(N+l) 



,N+l,N+l,app (v N+l n \ , ,JV+l,JV+l,appr / (JV+l) 



|y^+l)|2_ 2 "i 



(30) 



From the induction hypothesis, we extract the error term of order N + 1 of Eq^ ,n ^ : 

MAO 

~ (AT) ^ ~ (N,n) 
n=N+l 

~ (N,n) 

with Eqj ' admissible of degree n. We then rewrite: 



L + (i? e (^ +1) )) + A|(6f +1 )). % .Q' 



L-(Im(Re{T<- N+1) )) + A,(mf +1) )AQ] + ^? 

fe 2 (iV) 

+jErr ( 1 ) + ^ rr ( 2 ) + £; rr (. 3 ) + £ rr ( 4 ) + Eq 

3 3 3 3 Z <■ 



~^(JV,JV+1) 



,N+l,N+l,app n 
3 |Q| 2 ^ 



(JV,n) 



(31) 



n=N+2 



Step 2. Computation of the order N + l corrections. 

Let us now compute the order N + 1 corrections. From (VQ,yQ) = — \\Q\\ 2 L2 we may 



chose b^ N+1 ^ such that 



N- 

\Q\' 



1>jv+1 ' orp q),vq) =0 
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From Lemma [21 we now can solve for a Re(T, 



L + (Re(Tj N+1) )) + (b) 



Vj)Q 



j 

+ Re 



satisfying: 



Eq 



(N,N+1) 
j 



,N+l,N+l,app 

\Q\ 2 



Q 



Similarily, from the subcritical relation (AQ,Q) = |||Q||^2> we chose rrij^ satisfying 



-\ jmj 



(JV+l) 



AQ + Im(Eq 



~ (N,N+1) 



J 



,N+l,N+l,app 
9 ir)|2 



Q),Q)=o 



and thus we can solve for: 



L-(Im(T- )) + {nij 



)AQ 



+ Im 



Eq 



(N,N+l,m) 
j 



,N+l,N+l,app n 
\Q\ 2 t * 



0. 



Observe now that <^|^p 1,Ar+1,app <5 is admissible of degree N+l and thus from the induction 
hypothesis and Lemma [21 we just count the homogenity in a and find that Tj N+1 ^ is 
admissible of degree N + 1 and nij N € <Sjv+i, bj N+1 ^ G Sjv+i- 
Step 3. Estimate of the error. 

We now estimate the error terms in the right side of (|3ip . Note that the last term in 
the RHS of (|3ip is by constrution admissible of order > N + 2. It thus remains to consider 



the (-Errj fc) )!< fc < 4 . 



Err^p: Consider ([2 
from deg(yj N+ ^ 



From the construction of Vj N+1 \ Mj N+1 ^ and Bj N+1 \ and 



Q) > 1, deg{Mf ) ) > 2,deg{B i f ) ) > 2 and deg(T^ N+1) ] 
deg{mf +i> - - ^ {N+1) ^ kt , 1 cwC 1 ) 



), 



> iV + 1, EVrj is admissible of degree > iV + 2. 



Err 1 ?': Consider ([28]). We observe that given tp admissible of degree n, ^ are 
admissible of degree n while is admissible of degree n + l. We then use the 



8V (N+1) (N) 

fact that deg( — -kr — ) > 3 and the key degeneracy deg(B- ) > 2 to conclude that 



(2) 

Err) is admissible of degree > N + 2. 



Errf ] : Consider ([29]). Then deg(\V. 



together with Lemma [T] imply that Err - is admissible of degree > N + 2 



(AO 1 2 
(3) 



Q 2 ) > 1 and dec/(V, ( 



(JV+l) 



Q) > i 



ErA 4 ^: Consider (|30p . The key is to recall from (iii) of Lemma [T] that 



adds at least one to the degree of if). Using again deg(Vj 
that deg{Errf ) ) > N + 2. 



(JV+l) 



,N,app 



always 



Q) > 1, we conclude 



This concludes the induction step and the proof of Proposition [TJ 
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2.6 Definition of the approximate solution 



Now we explicitely describe how to extract a high order approximate solution to (pQ) from 
the profiles defined in Proposition [TJ 

Proposition 2 (Definition and accuracy of the approximation) Given N > 1 and 



parameters P = (a\, a 2 , fix, @2, Ai, A2) € 



5* \2 



(^ (Ar) )p(%) = V^\a 2 - aufo - 1 ,XiA2,V j ), 3 = 1,2 

be the object constructed from Proposition U\ 

Let t -» P(t) = (a 1 {t),a 2 (t),p l {t),(3 2 (t),X 1 (t),X 2 {t)) and 7 (i) = (71 (t), 72 (t)) 6e C 1 
functions on a time interval I = [to,tx], t\ < +00. Let a(t) = a 2 {t) — a\{t) and (3(t) = 
~ 0i(t) and assume that the following controls hold: Vi € I = [to,t\], 



i<MoI<|, Q 



< 



2 

A,- (to) 



<2||/3(to)| 
< Xj(t) < 2A i (i ), j = l,2. 



Define the approximate two soliton solution by: 

R( N \t,x) = Rf ( l, t Jx) = R[ N \t,x)+4 N) (t,x) 



(32) 
(33) 

(34) 



Xj(t) 



Then R( N ' satisfies 



i=i i 



S 



(N) 



A,(t) 



(35) 

where given by \18\) encodes the dynamical system of the geometrical parameters and 



error terms and satisfies: 
Vx£R N , \^ {N \t,x)\<- 



N+l 



-C f JV ) A 1 ,A2,^ 1 ,^ 2 (ll a! - Q! i(*)ll+ll a; -«a(*)ll) 



(36) 



Proof. According to the computations leading to the derivation of the approximate 
equation (117j) . there holds the explicit expression of 



(iV) 
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(37) 



2 1 



.Eg 



(JV) 



Aj (t) 



k N,app 



y W|2 



\H-i 1 \ A, 



A,- , a 

—vM-iy— 



A 



•i+i 



(JV) 



.T- 



-Oj(t) N 



A,(t) 
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with Eqf' given by (122]) . 

The first term in the RHS of (|37p is estimated from the exponential localization in 
space of the profile Rj ■ 



~ C|l ^e(< ) < y ) ll ^ (|jR ^ )| + l ^ 

< lli? (7V) i? (Ar) ll 3 e _C7jvr ' A i .^a.Pi.Pa CII*-"! (*) ll+U^-ota (t) ||) 
— II 1 2 

< g-Civ.ALAa.^^dlaWID+lk-aiWII + lk'-^WII 



W| 



where we used the standard Hardy-Littlewood-Sobolev inequalities. 

According to Proposition [H the term Eq^ in the RHS of (I37j) is admissible of degree 
> (iV + 1) and estimate ()36p follows for this term. 

It remains to estimate the third term in the RHS of (|37p that is the error induced by 
the dipole expansion of the gravitational field created by vj^-y . We claim that: 

2 

V: 



,N,app 

\V (N) \ 
v j+i I 



A; 



'3+1 



AN) , 2 

I \ Aj + l 



< CnmmaaO- + \\Vj\ 



\N+1 



A 

^•fell< i„ 



' v; + (-i) J 



Aj+i 



1, „ „^||a|| 



X \\N+1 + 1 A J || W ||>^ 



(38) 



The exponential localization of vj in yj then implies 

2 



k N,app 



A, 



V 3 + l V A j + l 



A 



V {N) < - 

3 — 



-CV.A-^Aa, ,3i,/32 fell 



a 



\N+1 



which concludes the proof of (|36|) . 



Proof of (|38p : We first observe from exponential localization of and the defini- 

tion ([To]) of the dipole expansion that: Vyj E 



,N,app ( ■ 



+ 



■3+1 



\ v j+i\ 2 V A 



A 



■i+i 



< Cn^maaO- + \\yj\ 



\N+1 



We may thus restrict our attention to the region Aj||j/j-|| < We then split: 

2 



A, 



A,- , a 
1 -Vj + (-1) 



Aj+i 



a 2 r W^(0\ 2 



4irX j+1 J R 3 \\a - (-l)i(A i+ i^ - A^j) 
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4vrA 



■i+i 



l^(OI 2 



3 . +1 ||4||<M \\a-(-iy(\ j+1 Z-\ jyj ) 



+ 



\a-(-iy(X J+1 (-X jyj )\ 

(N) 



From the exponential localization in space of V^ +1 and Hardy-Littlewood-Sobolev: 



l*#?(Ol a 



l?ll> J 



\a-(-iy(\ j+1 c-\ jyj ) 



< C 



Ai,A2 



V {N) 1 



2 

8 

LI 



Now for AjllyjH < and Aj+i||£|| < we use (I14p that is the uniform estimate: 



VH>i, V||C||< 



JV 



l«-CII 



fc=i 



to derive: 



4^A 



i+l 



i^iS (or 



(i + ||CH)" +1 



a 



IJV+1 



J + J.UMI — 10 

'1 + ll„,lhJV+l 



,Af,app / x 



IJV+1 



<i ) (Oi 2 (i + iieiir +i ^ 



AT+l, 



and (|3~8|) from the exponential localization of vj+± . This concludes the proof of (f3~8j) . 
This concludes the proof of Proposition [2j 



3 Reduction of the problem 

In this section, we reduce the proof of Theorem [1] to a bootstrap argument i.e. Lemma 
The reduction proceeds into two steps: a compactness argument assuming a uniform 
estimate (section 13. ip and a proof of the uniform estimate assuming that the bootstrap 
holds (section f3.2|) . 

This approach follows the strategy developped in [23], [E] and allows us to use an 
energy type method for the construction. 

Let (Af°, AEf) e (R* + ) 2 and a solution 

P°°(t) = (a?(t),a?(t),P<?(t),f3?(t),\?,\?) 

of the Newtonian two body problem: 

f af = 2/3f, a? = 213?, a°° = af - af t 

\ ooo _ Mi / ooo__Mi ^ 
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satisfying 

\?a?(t) + \?a?(t) = 0, (a i )^ ) (t) = 0, (40) 
and either of hyperbolic type with any masses 



4" " 4?r VAf° A? / ||a c 



or of parabolic type with the restriction of equal masses: 

E = and \™ = \™ = A°°. 

Recall that the two particles move away from each other with a speed depending on 
the regime 

Hyperbolic regime: \\cP°(t)\\ ~ C P °ot, ||/3°°(t)|| ~ C> 0, (41) 

+oo +oo 

Parabolic regime: ||a°°(t)|| ~ C P ood, \\P°°(t)\\ ~ — . (42) 

+oo +00 £3 

By rotation invariance, in addition to the assumption that the P°° trajectory takes place 
in the plane {^(3) = 0}, we may also require without loss of generality that in this plane, 
the line {a>(i) = 0} strictly separates the two particles, in the following sense 

liminf ■ ei > 0, liminf . (_£) > 0, (43) 

where \i = 1 or |, depending on the regime. 

3.1 A compactness argument 

Let us first show that the system of modulation equations inherited from Proposition [1] 
and (|18|) may be solved from infinity with asymptotic behavior given by P°°. 



Lemma 3 (Integration from infinity of the refined modulation equations) 

LetMf\ B p be given by 1121]) in Proposition^ There exists a solution 



of the system 

< a[ N) = 2(i[ N \ 4 N) =2& N \ a W=aW-af\ fjW = /?f > - /?f > , (44) 
1 ^ = <V VU^Af ) ), = ^(aW.^.f.f), 

on [Tg, +00) /or Tq large enough, satisfying in addition: 
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1. Hyperbolic regime: 



|i^(f)_ P W(f)||<Sfl. 

t2 



2. Parabolic regime: Assume AJ° = X% . Then, Vj = 1,2, 



fc=i 



K°)(fc)(*) 1 



+ !(«}%)(*)! + (t) - + |A^(t) - A 



< 



l ) 

is 



Proof of Lemma [3J, See Appendix A. 
Remark 3 In particular, there holds using \^ fi42\ ) as t — > +oo; 

* {A °(i)|| ~ Cpcot, ||/3 (iV) (t)|| ~ Cpoo > in i/ie hyperbolic 



\ot 



case, 



\a^>{t)\\ ~ CW S , 11/3^(4)11 + ||/3r j (4)|| < in the parabolic case. 

+°° £3 



(45) 
(46) 



We claim the following uniform backwards estimate, which is the heart of the proof of 
Theorem [TJ 

Proposition 3 (Uniform estimates) For N > 1 large enough the following holds. Let 
p( N ^(t) be the solution obtained from Lemma\^and let ^ N \t) = (7^ (t), 73 (*)) satisfy 



J = 1,2, 7 f } 



(JV),,2 fJV) 



(A w )2 



Lei i? 



(TV) 



Lei a sequence T n — > +00 and u ra (t) 6e £/ie solution to 

j id t U n + Au n - (p\ Un \2U n = 0, (/>| Un |2 = ~4^pn * |«n| 2 , 
\ n n (T n ) = ^ ) (T n ). 



(47) 



T/ien i/iere exists Tq = Tq(N) independent of n such that: 
Vn > 1, Vi € [T ,T n ], u n (t)-R 



(AT) 

PW(t) 



< 



ff 1 



(48) 
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Let us prove Theorem [I] as a direct consequence of Proposition [3j 

Proof of Theorem [1] This proof follows the lines of the one of Theorem 1 in [21] . 
By time translation invariance, we suppose without loss of generality that To = in 
Proposition [3j 

By (HH|), it is clear that there exists C > such that for any n, t G [0, T n ], 

\\u n (t)\\ m <C. (49) 

Now, we claim the following compactness result. 

Lemma 4 There exists Uq € H 1 ^ 3 ) and a subsequence {u n >) of (u n ) such that 

Un'(0) — > Uq in L 2 (R 3 ) as n — > +oo. 

We sketch the proof of Lemma [U Note that local convergence in I? for a subsequence 
is obvious by the H 1 bounds. The lemma thus follows from the following uniform L 2 
localization property 

Ve > 0, 3A = A (e ), such that Vn > 1, / \u n (0,x)\ 2 dx < e. (50) 

J\\x\\>A 

Note finally that (|50|) is a consequence of (j48j) (here To = 0) and simple computations. 
See the proof of Lemma 2 in [21 J for more details. 

Consider now the global H 1 solution U(t) of 

iUt + AU — 4>\u\iU = 0, ^ = -3-^*1^, 

(t,i)eixi 3 , [/(o) = c/ . 

Fix t > 0. For n large enough, we have T n > t and by continuous dependence of the 
solution of ([1]) on the initial data in L 2 (M 3 ) (see [1]), we have 



u n ,(t) -> U(t) in L 2 (R 6 ) asm +oo. 
Since u n >(t) — R^ N) ^ is uniformly bounded in H x , it follows that 

u n ,(t) - RW - U(t) - in H^R 3 ) as n ^ +oo, 



and thus, by (|48j) . we obtain 

^ > n \\TT(A - R (7V) )(t 



Vi>0, (52) 

W t 8 

Therefore, Theorem [T] follows from the definition of itf^,, in Proposition [2] and the 



properties of P^ N \t) described in Lemma El 
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3.2 The bootstrap argument 

We now focus onto the proof of Proposition [3] which relies on a bootstrap argument based 
on energy type estimates. 

Proof of Proposition [3l Let N large enough to be chosen and let R^r)u\ j(n)m De 
chosen as in the statement of Proposition [3j 

Let u n (t) be the solution to (f47l) . As long as the flow evolves close to the two solitary 
wave flow, we may from standard arguments (see e.g. [19J, Lemma 3) use the implicit 
function theorem to prove the existence of a unique C 1 geometrical decomposition of the 
flow: 

u n (t,x) = Rpty^ix) +e(t,x) 

(see AMI) with 

P(t) = (a 1 (t),a 2 (t),f3 1 (t),p 2 (t),\ 1 (t),\ 2 (t)), 7 (t) = ( 7l (i), 72 (i)) 
chosen such that the following orthogonality conditions hold: for j = 1,2, let 

Bj(t,yj) = Xj(t)e(t, \j(t) yj + a J .(t)) e ^(*)-^(*)-(^(*)%+^( t )), (53) 

then: 

Re( £j (t),Q) = Re{sj{t),^~Q) = Im(sj(t),AQ) = Im{£j(t),VQ) = 0. (54) 
Observe from (|47f) that: 

P(T n )=pW(T n ), 7 (T n )= 7 (A °(T n ), e(T n ) = 0. (55) 

We claim the following bootstrap lemma which is the core of the argument: 

Lemma 5 (Bootstrap) For N large enough, there exists Tq(N) > large enough such 
for allT* e[T ,T n ], if: 

Vt G [T* , T n ] 

then: 

Vi € [T* , T n ] , 



{ \\e{t)\\ H i < -1, 

>2 



E 
E 



2 



|Ai(t)-Af5(t)H-||/%-^| 



< 



< 



JV ! 



(56) 



Ej=i 
Ej=i 



JV ) 

2t^F 



\\ 3 (t)-\f\t)\ + w 3 -p) 

\^t)-lf\t)\ + \\a 3 -a^\ 



< 



2? 1 



JV 5 



< 



v 

2i"B" 



(57) 



25 



From a standard continuity argument, Lemma [5] proves that (|57p is indeed satisfied on 
[T ,T n ]. Thus, Vi £ [T ,T„], 

IU, m - R (7V) IKL W-bW II , II F?(iV) „(iV) || _C 

H«nW ^pW(t), 7 W(t) II — II "iW ^(t)^*) II H ri P(t),7(t) ^pw^iv)^)!! ^ iv 

which implies (|48|) on [To , T n ] and concludes the proof of Proposition 
The rest of the paper is devoted to the proof of Lemma [5j 

4 Control of the dynamics in the bootstrap regime 

This section is devoted to the proof of Lemma [5] which relies on energy type - and not 
dispersive - estimates on e. 

4.1 Control of the modulation parameters 

We focus onto the proof of the control of the modulation parameters (|57p which is a 
consequence of our specific choice of orthogonality conditions (|54p . 



We follow the notation of the proof of Proposition [3l We simplify the notation R{t) = 
PH)f(t) = + ^(t)- First observe that (|4"5j) . ([16]) and the bootstrap assumption 

(|56|) imply on [To, TV] for To large enough: 

||a(i)|| ~ Cp°ot, \\P(t)\\ ~ Cpoo >0 in the hyperbolic case, (58) 

+oo +oo 

2 C' 

\\a(t)\\ ~ Cpoctt, (t) || + ||/? 2 (t) || < —t 21 in the parabolic case. (59) 
Step 1. Structure of the modulation equations. 

Let us denote by Mod(t) the errors of the modulation equations with respect to the 
law predicted by the dynamical system P^: 

2 1 

Mod{t) = ||dj - 2^|| + |! Xj - Mf> || + ||4- - Bf> || + \jj + - 11^- 1| 2 - $j ■ a, | (60) 

3=1 J 

where iwf } = <V) = ilff>(a, AA 1; A 2 ), = = B f\a, (3, \ u \ 2 ). 

Now, we write the equation for e. From (|35p and e = u n — R, there holds: 

^ t £ + A£-2^ e{el) i?-^| i?2 |e = *W+jV(e) (61) 



V- 1 \ q (N)] f X-ajjt) ^ 
^ ( A,-(t) )e 



i=i A ^ 



-i"fj(t)+i0j{i)-x 
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with 

M{e) = <j)\ A 2R + 2<t) Rem e + 0| E pe. (62) 

Next, we compute the modulation equations using the orthogonality conditions (|54p . It 
suffices to take the inner product of (|61|) by four given functions exponentially localized 
around each center of mass (ay(i))j=i,2- Let Qj = Qj(x) be a given function with some 
decay at infinity and let 



G(t,x) = 

X 3 



1 / x-ajit) 



A, 



We compute from (|6ip after an integration by parts: 

Jm y eel =Re j e \id t Q + A© - 2^ e(e7!) i? - 0j^jG 



eft 



-Jte( | + AT( £ )]9) + Y J Re([y 4 

j=i \ i 



S 



(TV) 



We now use the fact that the approximate modulation equations for scaling and Galilei 
are both of order 0( ipmr) and the localization of around ay to derive: 



f_L + Mod(i)j e -C(||x-ai||+[| a! -ca||) i 
where we defined the linear operator close to vj N ^ by: 



Lj9j = -AQj + Oj + 2<f> 



Re{&M N) ) 3 



J 



47rA, + i 1 1 or I 



Qj (63) 



Note that we have approximated the field created by the by the first order of its 

dipole expansion thanks to the space localization of Qj up to an O(jnjjp) term. We thus 
obtain the preliminary system of modulation equations: 



j t { I'm I ;B ^ R< I : 



A4L 3 



i\ 2 ;dt@j - LjQj 



{jJ-) e - i l3( t )+ i Pj{ t )- X 



(64) 



° N +¥^ + C N \\e» 2 



Ml^ 1 Hall 2 



n 1 



Step 2. Approximate null space. 
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If we replace Lj by the exact linear operator close to the ground state, then the choice of 
orthogonality conditions (j54[) corresponds to the null space of L and generates modulation 
equations that are quadratic in e. In our setting, we claim that the choice (1311) reproduces 
the null space structure for Lj up to an O(-rpp-) which is enough for our analysis. 



Indeed, let us rewrite the equation of (f22j) as follows: 



Vf°l 2 47rA i+ i||«| 



V (N) = O 



-C\\x— a,j 



(65) 



We then let the group of symmetry of (pQ) act on this equation and differentiate with 
respect to the symmetry parameters to get: 



-C\\x—a 



J ", (phase), 



(66) 



-1 + 



*3 



47rAj + i \\a 



V™ + O 



L J (VV} N) )=0 



a 



(■ji) e-OH*-^^ ( sca i ing ) ^ ( 67 ) 

(68) 



-2V^ (7V) + O 



CIN oLj\\^ {translation 
1 



a: 



(Galilei). 



(69) 



We now compute the modulation equations from (|54p by using (|64p with successively 

. Observ* 

(Mod(t) 



<dj = ivj N \ AVj , wj N \ iyjVj . Observe that for these four functions, 



id t Qj = O 



+ 



1 



a 



a 



e -c\\x- aj \ 



and thus (j66j) . (|67j) . (|68 |) . (|69 |) yield: 

2 



Mod(i) <J2 Re f 

3=1 



+ 



CMod(t) C_ 



+ 



N 



a 



\N+l 



+ 



IhIh 1 



a 



+ C N \\e\\ 2 Hl . (70) 



We now use the dipole expansion on the term <ft Re (Q-R-} > as i n section 12.21 the orthogo- 
nality condition Re(J eRj+i) = from (|54p and the nondegeneracy relations J AQQ ^ 0, 
J yiQdiQ ^ to conclude: 



Re I e ^ J Re(e7?7) i? j+l 



Hall 



Re( J eR j+1 )+Q 



O 



\ e \\m 



Injecting this into (|70p yields 



Mod(t) < 



Cn 



a 



\N+1 



+ 



\ £ \\m 
||a|| 2 



+ CjvlMliyi- 



(71) 
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Step 3. Integration of the modulation equations. 

We now prove the estimates on the parameters of (I57p as a simple consequence of (171 fl . 
We inject (jBT)]) . ([55]) . ([55]) into dTTJ) to get for T* > T (iV) large enough, 



Vt G [T*,T n ], Mod(t) < 



(72) 



From (|2ip . the definition of A^- , (see Lemma [3]) and the explicit structure of 

M j N \ B j N) of de g ree > 2 > the re holds: 

lA.-Afl + ||ft-/f>|| 



< 4 + - M f , ( P<N, )I + » B f '( p ) - B f '( P<iV, )ll 

2 



1 C 



(AT) 113 



a — a 



+ 



c 



a 



(AT) 112 



(AT), 



fe=l 



< 



c 



9 , AT • 



Hence from ([55]) : Vt G [T* , T n ] , 



(AT) 



< 



c 



1+4 



(73) 



We now integrate the translation parameter and find from (|72p and ()73|) : 



| a j - a 



M"<C||#(t)-/5^'-' ' ' 

tT AT* 



+ — < 



1+4 



from which: Vt G [T*,T n ] 



\otj(t) — 1 



< 



C 



1 "' 7v 2 r« 

Next, we estimate the phase parameters from (1721 . (|73l) and ([74]) : 



(74) 



< C |^-Af | 



%-ft (A "ll+*l 



(TV) 112 



< 



c 



ivt i+ - 

and thus: Vt G [T*,T n ], 



l7i(*)- 7f } (t)l 



N 



(75) 



Finally, ([73]) . ([71}) and ([75]) now yield the estimate for the parameters in ([57]) for N large 
enough. 
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4.2 Control of e in H 1 

We now focus onto the proof of the bootstrap estimate (I57p on e. 
Let us consider the energy conservation: 



2H = 2H(R + e) = 2H(R) -2Re(e,AR- (j>\ Rl 2R) 

+ J |Ve| 2 + J <f> lR{2 \e\ 2 - 2 J \V<p Rem \ 2 + 2 j <f> Re[£R) \e\ 2 -\ J lW k , 2 | 2 . 

We add to the quadratic and higher order terms of this conservation law a localization 
of the two other conservation laws, namely I? norm and momentum conservation, on 
each solitary wave. We obtain a functional which is almost conserved by the linear flow 
governing the equation for e (|6ip (see [19] and [21] for similar technique). 

Let ^l{x) = (recall the notation x = (xh\, Xt 2 ), ^(3))) denote a smooth non- 

negative cutoff function such that ^(x) = 1 if xm < —1 and ty(x) = if arm > 1. We 
distinguish the two cases of the asympotic dynamics P°° being a hyperbola or a parabola: 

1. Hyperbolic regime: According to (|58p . we may find a constant C = Cpoa such that 

ipi{t, x) = Mt, x) = l- x) 

localizes around each center of mass otj(t), that is ipj{x) = 1 around ay(t) (i.e. in 
some region of the form \xm — (ay)(l)(^)l < Ci)- We then let 

Cj(t,x) = ij)j{t,x). 

2. Parabolic regime: According to (|59|) , we may find a constant C = Cpoo such that 

iJ> 1 (t,x) = V(-A r ), ^ 2 (t,x) = l-'0i(t,x) 

2 

localizes around ay(i) (|aim — (ay)n)(i)| < Ct%). In this case, we let 

1 



Cj(t,x) 



2 



The point is that the common notation Q will allow us to perform part of the proof for 
the two regimes at once. 

In both cases, we define a functional in e(t) which aims at localizing in space the global 
conservation laws: 

|2 of \r-rj, 12,0(1 |_|2 1 



(76) 



G(e(t)) = \Ve\'+ /0 |Rp |^-2 / |V^ e( ^| 2 + 2 /^ e( ^|^-- / |V0, 




j y Cjkl 2 - 2/3j J tpj ■ Im(Vee) 
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Note that for the parabolic case, the L term does not contain a cut-off (the proof would 
not work, this is where we will use \ f = ASf 5 ), while the momemtum term does contain a 
cut-off. For the hyperbolic case, both term contain a cut-off so that any A°° is possible. 

We first claim the following coercivity property of the linearized energy which is a 
standard consequence of the variational characterization of the ground state and the choice 
of the orthogonality conditions (|54p , together with the nondegeneracy of the kernel of the 
linearized operator L close to Q as proved by Lenzmann [15J: 

Lemma 6 (Coercivity of the linearized energy) There exists a constant cq > such 
that under the orthogonality conditions \54\ ), for any t >Tq large enough, 

g(e(t))>c \\e(t)\\ 2 Hl . (77) 

Proof of Lemma [BJ. See Appendix B. 

We now claim the following energy estimates on e which is the core of the proof: 

Lemma 7 (Global energy estimate on e) There holds for N , T* > To large enough: 
Vie [T*,T n ], 

d 



c, W" 1 



with constant C > independent of N. 



Remark 4 Observe that (7?) and (73|) formally imply a structure: 

d 



dt" £ 



1 2 

iH 1 



. r( _^_ + MMIi). (79:i 



ii £ (*)ii 2 

A fundamental difficulty here is the gain j only in the control of the error terms 0{ j^ 2 -) 

mostly induced by the localization of the conservation laws. This size comes in the hyper- 
bolic case from the distance \\a{t)\\ ~ t. In the parabolic case, the situation is much worse 

2 

||a(t)|| ~ ts and hence we need to restrict to the symmetric case, see estimates &83\) . t90\) . 
\91\) below. Reintegrating |7ff|) backwards from infinity, we see that we may bootstrap an in- 



formation e \\ m < -^r only if N is large enough, hence the requirement of an approximate 
solution to a large enough order. 

Let us postpone the proof of Lemma [7] and complete the proof of Lemma \E\ 
Proof of Lemma O 

We may now conclude the proof of (|57p . Here the key is the fact that the constant N 
in (|78p may be fixed as large as we want. We integrate (I78p in time so that from ([55]) and 
dS6D: Vt G [T*,T n ], 

\Q(e{t))\ < / -tttt + — — Mr < -rr + —w < 



T T ^ 2 / t Nt 2 Nt 2 
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for t > T* > Tq(N) large enough, where C is independent of N. From (|77p and (|57p . this 
implies: 



colkWII^ < 



and (I57p follows by fixing N such that 

Nco > 8C. 

This concludes the proof of the bootstrap Lemma [5] assuming Lemma [71 

|| E ||2 

We now turn to the proof of Lemma [7J We point out the key point of gaining the — j 2 — 
degeneracy in the RHS of ([75]) . the power t being sharply enough to close the bootstrap 
estimate as in the proof of Lemma [5j 

Proof of Lemma [7J. Let us split: 



Q = Qi + Q2 + Q3 



with 



Qi = J |ve| 2 + J <WN 2 " 2 / I v <Wr)I 2 + 2 / ' 



'Re(eR) I 



|V4 



Step 1. Estimate on Q\, We compute: 



d „ , . 



-2im / i^e 



+ 2 / |e 



Ae - 0| fl |ae - 2</> ite(e7 j ) # - 2<^ e(e7?) £ - |e | 2 i? - 0| £ pe 

S) + 4jRe / ^t^ReieR) + 2Re / ^0| £ P- 



J Re(dtRR) 

We now observe from (|34|) and from computations similar to (|1(J|) that: Vj = 1,2 
1 



d t Rj(t,x) 



AU 



and thus from J56J, and (gTj, (02}: Vj = 1, 2, 



-2/3,-V^ + z^ + ll^ll 2 ^ 



-Cjv||a;-aj 



(80) 
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We now use the localization properties of Rj around ay, inject the equation for e (|6ip and 
get using (f80l) . the estimate of vl/W ([36jl . the estimate on the geometrical parameters (|72l) 
and the bootstrap estimate (|56p for cubic and higher terms in e to derive: 



d r ^ r / ggljg(^ , \\e(t)\\ 2 Hl \ 

7 = 1 



with 



'/,V.,',V />'./>>,: 



(81) 



Step 2. Estimate on C?2- We claim that: 



d_ 

dt 



+ O 



(82) 



We argue differently in the hyperbolic and parabolic regimes: 

1. Hyperbolic case: From (|58p . the distance beween the two solitary waves in the hyper- 
bolic regime is lower bounded: ||a(i)|| > Ct, C > 0, which allows the estimates: 



NCj(t)\\ + \dtCj\ < 



C 



(83) 



Moreover, + |Aj(t)| < ^ from the equations of motion. Thus, from (|6ip using also 



A 2 



+ lift 



ON 



^ + ll/?# 



Im(id t e(j£) + O 



2 / Jm(VCj • Vee) + 4 / Im(eQR)(f> 



+ (^^ + M|^ + l|£(t)lliflM ^ ) ), 



(84) 



and ([82]) then follows from (fflj) . 
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C 

t3 



2. Parabolic case: From ||a(i)|| ~ is in this case, and thus + |Aj(i)| < 

Since ( = \ , using (|TTj) and arguing as above, we get 

^ _ 2 £ ( i + lift ||») / M.^^ + o + 

+ O ([|A a - A 2 | + - ||A|||] 11^)11^) + O ( ^f^ + 
We now claim that: 

|Ai-A 2 | + |||/3 1 ||-||/3 2 ||| <^ (85) 

and ([82]) follows. 



Proof of (185]): We estimate from {60]), ([72} and Xj(T n ) = X ( f ] (T n ) from (g7|: Vi £ 
Po,T n ], 

|Ai(t)-A 2 (t)| < |A 1 (t)-Ai 7V) (t)| + |Ai 7V) (t)-Af ) (i)| + |A 2 (t)-Af ) (t)| 
< / Mod(r)dr + |A^ ) (t)-Af ) (t)| 



< ^ + \x[ N \t)-Xi N \t)\. (86) 



We now estimate from the definition of Pat and (1251) : 



< 



\W (N) \\ u(iV)_ A (iV) |+ 1 



laW 1 2 1 llaW 



(87) 



Now observe from Lemma [3] and the symmetry assumption on the asymptotic trajectory 
in the parabolic case that 

lira xf\t) = X? = X?= lim Af>(i) 

t— >+oo t^+oo 

and thus the time integration of ([87]) with ||a^(i)[| ~ i 1 , \\P^ N ^\\ ~ i^a yield a first 
estimate |Aj (t) — A 2 ^(i)| < 4r- Reinjecting this bound into ([87]) we eventually derive: 
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lA^-A^)!^ 



which together with ([86]) yields ([85]) for Ai — A 2 |. 
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We argue similarily for (3. From fl60|), £73]) and 0j(T n ) = ^ N) (T n ), there holds: Vt G 

+/? 2 (t)ii < wm - p[ N \t)\\ + \\p[ N \t)+pi N \t)\\ + n^(t) -pi N \t)\\ 



< 



i + nM Ar) (t)+/3f ) 



We now estimate from the definition of Pm and 



|A 



(AT) .(JV) 

1 — ^9 



(JV) 112 



+ 



a: 



(TV) 113 



~ t 2 



where we used (|88j) in the last step. Integrating this with lim^ +00 (3\ 
from Lemma [3] yields: 



(/V) 



limt- 



(TV) 



~ t 



i N) +pi N) \ 



which together with (|89j) concludes the proof of (|85|) . 



Step 3. Estimate on £3. First we have, since 



j t j \-2Im{^ ■ Vet)] 



4 / Re[ id t e 



< Ct~3. 



J 2Im(f3jd t ^j ■ Vee) + O 



Ht)\\h 



We now observe in the hyperbolic case that (|58|) implies: 

C 



f(t)||||V^(t)||L<» < 



t 



(90) 



A fundamental observation is that the same estimate holds in the parabolic case. Indeed, 
from fl59]): 

C „jjv),.,„ . C 



and hence 



< — + w j , m<- 

t 8 t3 



(91) 

This means that the localization in space of the kinetic momentum is harmless in the 
parabolic case even though both solitons are not far apart thanks to the explicit decay in 
time of the Galilean parameter. Note that by the same argument 



c 
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From the previous observations and ([B"T]) . ([56]) . (f72|) . we derive 
d 



(i/ , -2Im{p,il>j ■ Vee)} 



-V-(^)£ + ^Vi-Ve 



-4 y Pe I «d t e 

-4 /" Pe ([2<^ e(e7!) P + • Ve) + O 



+ 



H 1 



t 



( C N \\e(t)\\ m \\e(t)\\ 



2 

m 



and thus from the localization properties of ipj and Rj: 

jGz{t) = ^8 ( Re{el3 r VB~)(t) Re{eli) + ± ! klV^.v^Hj) (92) 
j=i ^ ■> 

+ / gjy||eft)l|gi , I|e(*)llgi 



Step 4. Conclusion. We add up (HH), ([82]) and |)2} to get exactly: 



d rM n ( C N \\e(t)\\ H i Hey* /!!/,-! 
and ([78]) follows. 

This concludes the proof of Lemma [7J 

A Proof of Lemma [3] 
A.l Hyperbolic Dynamics 

Let P°°(t) be chosen as at the beginning of section 3 in the hyperbolic regime. Let 
< e < g. Let To > 1 large enough to be defined later. We define fl e ,T = ^ the set 
of functions P(t) = (a\(t), ct2(t), Pi(t), fait), M(t), ^2(t)) defined on [To,+oo) such that 
|||P||| < 1 where the norm |||.||| is defined as follows 

2 

|||Pj|| = ^ SU P f^'IM*) - +* 2 " f WPj® - /fif (*)ll + |Ai(t) - Af|| • 

i=1 te[T ,oo) J 

Now, we consider the map V which corresponds to solve the dynamical system 
where the functions Bj N \ Mj N ^ are defined in Proposition [TJ For P(t) € ft, we set 

FP(t) = (r ai (t),ra 2 (t),r 0^)^02(1)^X^)^X2(1)), *E [To,oo), 
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where (as usual a(t) = 0:2(2) — oi(t), (3{t) = P 2 (t) — (3\{t)) 

/oo 
2(-(3 J (r)+f3]°(T))dT, 

/oo 

roo 

TX^t) = Af - J Mf\a(r),(3(r),X 1 (r),X 2 (r))dT. 

Now we claim 

Lemma 8 There exists Tq large enough such that T maps into f2. Moreover, we have 
the contraction property: let P a , P& £ £1, then 

|||rp tt -rn||| < i|||p a -p b |||. 

We note that this lemma immediately implies Lemma [3] via a standard Banach fixed 
point argument. 

Proof. We show that V maps VL into f2. The contraction property follows in identical 
manner by forming the difference equations. 

First, we consider Taj. Here we have, by |||P||| < 1, 

/■OO 

\\V aj (t) - a?(t)\\ = II J 2(/3,(r) - /3f (r))dr\\ < t^+i. 

Thus, choosing To large enough, we get 

||raj-(t)-af (t)\\ <t- 1+e . 
Next, we consider Tj3j{t). Note from (|39p and f|26|) that we have 

■,(2)/ 



P?>(t) = lim 0f(t) - / bf (a°°(T),/?*>(T),\f>,\%>)dT. 

t^oo J t J 

Thus, 

TPj{t)-Pf{t) = j^[bf\a^(T),^(T),Xr,Xf)-bf\a(T)^(T),X 1 (T),X 2 (r))]dT 
~ En= 3 JT 6f(a(T), /3(r), Ai (r), A 2 (r))dr. 

Using the expression of fr^ 2 ' 1 in (|26p . the estimates on b^ for n > 3 and the assumption 
|||P(£)!II < 1, we get (for e > small enough) 



"OO 



[bf(a™(T),(3™(T)A?A?)-bf(a(T),f3(T),Xi(T),X2(r))]dT 



< r 
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2 / 6f ) (a(r),/3(r),A 1 (r),A 2 (r),)dr 

n=3 ^ 



<t- 2 , 



whence we have shown \\Tf3j{t) - 0*°(t)\\ < r 2+ l. Choosing T large enough, for t > Tq, 
we infer 

||r&(t)-/3f(t)|| <t- 2+ i 

Finally, we consider TXj(t). Since rrJp € 52, we have \mp\ < i~ 2 , thus we get 
|rAj(t) - Xf\ < t~ 1+ i whence 



\rxj(t)-x™(t)\<t 

if To is large enough. This concludes the proof. 

A. 2 Parabolic Dynamics 

Here we prove the following: 



-l+e 



Lemma 9 (Existence of almost parabolic trajectories) Let A^ = A£° = A°° and 

P°°(i)(a<f (t),a% > (t),P% (t)iP2°( t )) be a solution of the two-body problem of parabolic 
type satisfying (Jty - There exists a solution (ax(t), ct2(t), (3i(t), ^(t), Ai(£), X2(t)) of (44\ ) 
defined for t > Tq with Tq large enough satisfying: \/j = 1,2, 



fc=i 



(«f )(*)(*) 



o 



0(fc) (*) 



+ |(aj)( 3 ) (01 + " /5f(*)ll + |A,-(t) - A°°| < 



C pa 
f3 



and 



\^i(t)a 2 (t) + X 2 (t)a 1 (t)\ -> as i -> +oc. 
In particular, the asymptotic trajectory is given by the same parabola like for P° 



(93) 
(94) 



Proof. Step 1. Leading order dynamics. 

First, we claim that the equation of motion of (i.e. (|44[) ) can be rewritten as 

follows: 



A; 



a ; 



+ 

2 ft, L 



IHr 



, J = 1,2 



w 



, c\a i v^2 

INF J =1 



2 C2a(Aj-A°°) 



IMI : 



(95) 



where co,ci,C2, C3 are explicit functions of (A°°,a,/3). Note that the same constants 
co,ci,C2,C3 appear for both j = 1,2, because of the assumption A 1 
over, the functions fj,gj satisfy the bound 



Vf = A 



§° = A°°. More- 



l/il<ll« 
\9j\ ^ II" 



I" 1 + \p\ + £ 2 =1 [|A, - A°°| + ||a||(Aj - A 00 ) 2 ] , 
I" 1 + S 2 =1 [|Aj - A°°| + \\a\KXj - A 00 ) 2 ] 



(96) 
(97) 
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as long as ||q|| > 1, \Xj - A°°| + \\/3\\ < 1. 

Indeed, the structure of (f95|) is deduced from a expansion of the gravitational field (flBT) 
to the order ||a|| -3 , and from structural symmetry properties: both Tj and Re(Tj )j = x t 2 
are chosen radially symmetric (see Remark [TJ. 

Let us introduce an intermediate system p( a PP)\ 



p(app) 



' ^(<w) = ^(app) = ^( app ) = c 3 (g( a ™\ f3^ ) 
1 2 1 1 (y.^ a P^ > ^ 11*^ ' 



a 



(app) 



~ 3 2tfT\ 0}' - 4 

and which implies the law for a^ app ^: 



■Japp) _ (-1)J 



(app) ClQ ("PP) 2C2o( app )( A(aPP) - A °°) 



CQQ- - - - I U1LE- ' ' - I 
\\ a (app) j|3 "T || a (app) p 



v(app) 113 



a 



cq a 



(app) 



a 



(app) 1 13 



+ 



C\OL 



(app) 



(app) 1 14 



+ 



2c 2 a^ pp \X (app ^ - A c 



a 



(app) 1 1 3 



(98) 



We claim that there exists a solution to p( app *) which asympotically converges to P°° in 
the sense of (1931) . 

Let us first observe from (j98[) that Ji(c^ app ) A d^ app -*) = and hence the movement is 
planar. We assume that cr app ) stays in the plane {^(3) = 0}. Choosing polar coordinates 
(r, 6) in this plane, a^ app \t) = (r(t) cos 9(t), r(t) sin 0(f), 0), we rewrite the system as 

' r 2 # = a > 0, 

\(app) _ £3r 



_co + <^ + 2c 2 (A( a Prt-A°°) 



We now use the classical change of variables u = - and rewrite the system as an ODE for 
u = u{6). From 



d dO d 



d 



dt dtde a ° u de' 



we compute: 



du 
- ao d9> 



d 2 u 



d\j 



de 2 



d9 



and hence the new system: 

( d\( a PP'> 
d0 



1 du 
'2 C3 d6»' 



1 «o(S# + «) = co - cm - 2c 2 (A^ - A°°). 
We are led to compare the solution a°° of 



oo\2/joo 



a 



(r°°) 

a oU(e°°) a ^ J 
u°°(o) = ^(o) = o. 



c i.e. r 00 ^ 00 ) 



c (l -cos(0°°)) 
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with a( app ) solution on (— #o,0] of the system: 

\ og(gf + u) = c - (ci - 2c 2 c 3 )n = c - C4U, n(0) = §(0) = 0. 

In fact, the system reduces to a linear one and the solution is explicit, but we do not need 
analytic formulas. The solution admits the following asymptotics near 6 = 0: 

u(6) = ^9 2 (l + 0(6)), r (6) = ^(1 + 0(8)), \^{9) = A°° - ^ 2 (1 + 0(9)). 



We compute the asymptotic time dependence using 

72. 

-0 

which yields: 



ao = r 2 e = ^-(l + 0{9)) 



0(t) = - (||) 3 (1 + 0(4)), r(t) = (^j 3 t§(l + 0(4;)), (99) 

X («pp)(t)-\™ = -c 3 (^f\ V§(i + o(l)). (100) 

Similarily, we compute: 

1 1 
= " (||) 3 (1 + 0(4)), r«(f) = (^) 5 tf (1 + 0(4)), (101) 

and hence the polar curves 

(r(t) cos(0(t)), r(t) sin(0(t))) and (r°°(t) cos(0°°(t)), r°°(i) sm(0°°(t))) 

are asympotic to each other as t — > +00. In particular, we have the asymptotic behavior 
of the cartesian coordinates of a^ app ^: 

< P) W=(^)"(l + 0(^)), agf ) (t) = -(^) 5 (1 + 0(1)), agf^O. 

(102) 

We now solve for (ai(t), a2(t)) solution to p( app ) by letting 

a a 

thanks to the symmetry in j = 1,2 of the system p( a PP\ and hence the a^ app \t) curve 
also has a center of mass fixed at the origin. Equations (f99~|) . (|100p . (|10ip now imply ([93]) 
for p( a PP\ 
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Step 2. Solving the full system. 

We now claim that we can find a solution to the full system such that 
2>$=i\\atj(t) - af pp) (t)\\ ^ as t -> +00. 
Recall that has been rewritten so that the following holds 



(104) 



; _ c 3 (q-/3) 3j(Q,/3,Ai,A 2 ) 

d = 2/3, 2/3 



, 3 



— n H 



1,2 

^2 



V 2 cgoCV-^!!) 1 /(«,AAi,A 2 ) • 
^ ' 7=1 up I [up > j 



1,2. 



where / = f 2 - fi satisfies ® . We let a(t) = a^ app \t) + S(t), \j(t) = X^ app \t) + fXj(t) 
and linearize the system around a^ app \ \( a PP\ Let us show how to bootstrap an estimate 



<4, ||«(t)||<4. 
|Mi(t)l + M*)| < j. 



(105) 



(106) 



for some large enough K; the claim then follows by a fixed point argument as in Appendix 
A.l, and thus will be omitted. 

The main linear term is expanded as follows 



a 



c> pp ) + 5 

\ a (app) + 



a. 



(app) 



\ a (app) 113 



+ 



1 



\ a {app) 113 



5-3 



a 



( a PP) . g 



\ a (app) ||2 



a 



(app) 



+ 



2 0(1) 



| a (app)||4 



We now observe from the explicit law (|102|) that 

a (app) . g 



\ a (app)\\2 



a 



(app) 



'(1) 



+ 



IH|o(i) 



2 

is 



and so 



Co a 

IMP 



(app) 



+ 



1 



45 



(app) 113 9^2 



2L I -u-^odi^ii-Hl^il 2 ) 



-25 



(3) 



i 2 +3 



Using the bootstrap (11Q5|) . (fMl) and (j 102[) . the linearized system takes the form: 



' Ai = 4-o(i), j = 1,2, 

is 

45 (1) 




-25, 
-25, 



(2) 
(3) 



with F(t) 



i 2 +3 



-0(1) 



(107) 
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An explicit solution is given by: 

ti f +oc F ^' iT> 

5(t) 



( 

V 



t3 



ts r°° F (2) (r)r3dr - ts J/ 00 F (2) (r)r3dr 
* f J T +0 ° ^(8)(r)ri(fr - t* / r +0 ° F (3) (r)r§dr 



from which f)105|) is obtained with better constant. 
To prove (|106j) . we have from p( N ^: 



< 



1 



(app) 112 



a 



a — a 



(app) I 



Uapp 



+ El A J- A(aPP) |)+°(l 



1 



a (app) 113 



) 



< ^s] =1 |A,-A^|) + io(l) 



which implies f)106|) by integration in time. 
We may now prove (|94h . Indeed, 



|A 2 (*)ai(t) + Ai(t)a 2 (*)| 



[A*2(*)II + Mi(*)||)t* + ||<5(t)|| < 



1 

"T 



where we used a 



(app) 



+ a 



(app) 



from (fT03D . 



This concludes the proof of Lemma [9l 



B Proof of Lemma [6] 

The proof follows exactly the same lines as the proof of Lemma 2.6 in [21]. We proceed 
into four steps. 

Step 1. First, we claim that there exists c > such that for any real- valued v E H 1 
(v,Q) = {v,yQ) = => (L + v,v)>c\\vf m , (108) 

(v,AQ) = (v,VQ) = => (L-v,v) > c\\v\\ 2 Hl . (109) 

Indeed, estimates f|108[) - (j!09j) follow from the arguments of [29], proof of Propositions 2.9 
and 2.10 and the fact that the kernels of L± are explicitely known from [15] (the kernel 
of L_ is spanned by Q and the kernel of is spanned by the components of VQ). 

Step 2. Second, we claim that for all e € H 1 satisfying the orthogonality conditions 
(f53|) - (f54"|) . and To large enough, we have 

Gj{e)>c \\e\\ 2 H1 , (110) 
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where 



9j(e) = j\Ve\ 2 + j^ Rj? \e\ 2 -2j\V<p Re{eR]) \ 2 

:| 2 - 2 J (i j ■ Im{Vee). (Ill) 



Let Qj(t,x) = jrfiQ C ^ ) e ~ il] {t)+ll3j {t) ' X ■ % the definition of vj N) in Proposition CQ 
the definition of Rj in Proposition [21 and the properties of the functions Tj we have 



IMI t3 



Thus, if we define 



J |Ve| 2 + J 0| QjP ki 2 -2 J \V<t> Re{eQj) \ 2 
|| 2 J J \e\ 2 -2 J (ij ■ Im{Vee), 



+ r 



(112) 



we have 



\Q^)-W\<^M\m- 

t3 



Then, by standard computations (see e.g. proof of Claim 7 in [2T]), we have 



Gj{e) = (L +J Re(e i ),i?e(e i )) + (L-Imfa^Imfe)) + 0(- 



1 



where are defined in §3§ . Thus, (fTTOj) follows from (fTU5j) - (fTUgj) and \\e\\ H i < C\\£j\\ H i. 
Step 3. The function t/jj being defined as in Section [3.2^ let 

12 



+ 



J 4>j\Ve\ 2 + J <^ |2 | £ | 2 -2 J \V4>Re { eR 3 )t 



Then, we claim that 



Gj, ps i(s) > c J i(}j(\Vs\ 



2 I |2\ ^ 11 ||2 

+ \ £ \ ) tIpIIj/i- 

t3 



(113) 



(114) 



Proof of (|114p . Let ij = e^Jtyj. Proceeding as in the proof of Claim 8 in [21] . we have 
/ - °\\4h < I ^|Ve| 2 < f + °\\e\\ 2 Hl 

J t3 J J 1 3 
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and ipjfij • ImiVee) = (3j ■ Im(Vejij] 
Next, we have 



and 



J 0|^i 2 kl 2 = y <P\R 3 \Atj\ 2 + j 0|^i a N 2 (i-V'i) 



+ 



(115) 
(116) 



Thus, by the decay properties of ip and Rj, we obtain 

J ^.| 2 |e| 2 -2 J \V^ Re{eRj) \ 2 - j <^.,2|e,f + 2 1 iV^-.^f 



<^IMI 2 

- 2 H £ lliJl- 
t3 



From these estimates, we obtain 



(117) 



By standard arguments (see e.g. [H]), since for t large ej- almost satisfies the same or- 
thogonality conditions as e, we obtain by step 2, Gj(£j) > cHe^H^i — TrlleUfj-i and (|114p 
follows. 

Step 4. Conclusion. We decompose G{e) as follows 

2 



i=i 



+ J 0|R|2_ ( | /?1 |2| + | ij2 |2)|e| 2 

- 2 j (|V0 ReM) | 2 - |V^ (£ ^ } | 2 - |V<^ e( ^ 2) | 2 ) 
+ 2 y ^ e(e «)|£| 2 -^ y |V0| £ | 2 | 2 . 



(118) 
(119) 
(120) 



By the decay properties of R\, R2, the term in (|118p is easily controlled by C/t. The terms 
in (|120p are nonlinear in e, and we obtain by the Hardy-Littlewood-Sobolev inequality 

y <pRe(eR)\e\ 2 + \f |V0| £ p| 2 < C{\\ef m + \\ef m ) < -w\\4m, 
where we have used (1561) for the last estimate. 
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Now, we estimate ()119p . First, by the decay properties of Rj, we have 

<i> M£mi) Re(e(t)R 2 (t)) <C [ e^-« 1 ( t )l[ £ (*,y)|e^l^WI|^ J z)| 1 -l- I dy ( te. 

j \y *i 

We cut this integral into two pieces depending on \y — z\\ > ^\\x\\ or ||y — z\\ < ^\\x\\. 
The first part, i.e. when \\y — z\\ > ^\\x\\, is easily estimated by ||e||^ 2 using Cauchy- 

Schwarz inequality. For the second part, we observe that if \\y — z\\ < ^||a|| then \\z— cks || > 
||«2 — "ill — ||y — Qti|| — \\z — y\\ > ^\\a\\ — \\y — ati\\ and so 

e --y\\y-ai(t)\\ e -"/\z-a 2 (t)\ < e --r\\y-ai(t)\\ e -^-y\\z-a 2 (t)\\ < g - |-y||y— cti (t)\\ & -\ \\a\\ _ 

Thus, 

e -7lll/-ai(t)ll | £ ( t) y )| e -7lk-«2(t)|| | £ ( t> z )| 1 ^ 

i/-*ll<£lMI \\y- z \\ 

< o e -ill«H f e-^-^WIlle^^He^^l^ < C 1 1 « 1 1 ^ e~ ^ 11 " 11 ||e||2 2 . 

J 2||<^ |]a|| 

Gathering these estimates, we obtain 

2 C 



Q{e) > ^Qj^ie) - —\\e\\ 2 m 



3=1 * S 

Using (I114p and f/'i + ip2 = 1, we find Q(e) > c\\e\\ 2 f r 1 — ^-\\e\\ 2 f r 1 which gives the desired 



result for Tq large enough. 
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